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Decoherence is the main process behind the quantum to classical transition. It is a purely quan- 
tum mechanical effect by which the system looses its ability to exhibit coherent behavior. The 
recent experimental observation of diffraction and interference patterns for large molecules raises 
some interesting questions. In this context, we identify possible agents of decoherence to take into 
account when modeling these experiments and study theirs visible (or not) effects on the interfer- 
ence pattern. Thereby, we present an analysis of matter wave interferometry in the presence of a 
dynamic quantum environment and study how much the visibility fringe is reduced and in which 
timescale the decoherence effects destroy the interference of massive objects. Finally, we apply our 
results to the experimental data reported on fuUerenes and cold neutrons. 
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I. MATTER WAVE INTERFEROMETRY AND DECOHERENCE 

Matter wave interferometers are based on quantum superpositions of spatially separated states of a single 
particle. However, as is well known, the concept of wave-particle duality is not applicable to a classical 
object because this kind of object never occupies macroscopically distinct states simultaneously. Then, by 
performing interference experiments with massive particles, in particular with those of increasing complexity, 
one can probe the borderline between these incompatible descriptions and shed some light on one of the corner 
stones of quantum physics. 

Matter wave interferometry has been largely studied in the last few years. Many theoretical studies have 
been done around the mesoscopic systems Mesoscopic objects are neither microscopic nor macroscopic. 

They are generally systems that can be described by a wavefunction, yet they are made up of a significant 
number of elementary constituents, such as atoms. Well-known examples these days are fullerene molecules 
Cgo and C70, which are expected to behave like classical particles. Nonetheless, the quantum interference 
of these molecules has been observed Q. In these experiments, also done with cold neutrons, thermally 
produced beams are collimated, diffracted by a grating, and then detected on a distant screen. The pattern 
so produced shows a typical interference profile of wave phenomena slightly attenuated. 

Usually, the main problem in the analysis of interference experiments is to establish exactly which are the 
causes for the loss of spatial coherence observed in the reduction of the visibility fringe of the interference 
pattern therein. Macroscopic quantum states are never isolated from their environments [Jj. They are not 
closed quantum systems, and therefore, they cannot behave according to the unitary quantum-mechanical 
rules. Consequently, these so often called "classical" systems suffer a loss of quantum coherence that is 
absorbed by the environment. This decoherence destroys quantum interferences. For our everyday world, 
the timescale at which the quantum interferences are destroyed is so small that, in the end, the observer is 
able to perceive only one outcome, i.e. a classical world. As far as we see, decoherence is the main process 
behind the quantum to classical transition. Formally, it is the dynamic suppression of the interference terms 
induced on subsystems due to the interaction with an environment. 
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In principle, some incoherence (lack of coherence) effects can be imputed to the passing of the particle 
through the slits, such as vibrations or Van der Waals interactions 0] or the difference in size of the slits 

0. In the present work we shall not consider such effects as they are in general negligible under suitable 
experimental conditions. We shall consider experiments where coherent states of massive particles are well 
prepared and the diameter of the particles arc smaller than the width of the slits in order to avoid the 
consideration of the above mentioned effects. 

In matter wave interferometry experiments, several losses of spatial coherence affect the particle beam 
during its evolution, with a consequence of a fringe visibility reduction of the detected intensity pattern. 
These dynamic decoherence effects can be imputed to collisions with the air molecules or thermal photons, 
for example. Formally, decoherence appears as soon as the partial waves (the wavefunction of the subsystem, 

1. e. massive particle) shift the environment into states orthogonal to each other. However, the loss of spatial 
quantum coherence can alternatively be explained by the effect of the environment over the partial waves, 
rather than how the waves affect the environment. It is a consequence of the entanglement between the 
system and its environment. The loss of spatial coherence can also be originated in the angular divergence, 
the non-monochromaticity of the beam and the randomness in the emission or arrival of the particle (mainly 
related to the experimental difficulty in the production of the same initial state for all the particles). This 
randomness gives raise to a fluctuating phase 4> ^md therefore, the interference term appears multiplied by a 
factor e"^. The effect can be directly related to the statistical character of (/>, in particular in situations where 
an external potential exerted on the partial waves is not static. We associate these effects to the dephasing 
process. Yet more importantly, any source of stochastic noise would create a decaying coefficient. In this 
way, the uncertainty in the phase produces a decaying term that tends to eliminate the interference pattern. 
This quantum suppression is due to the presence of a noisy environment coupled to the system and can be 
represented by the Feynman- Vernon influence functional formalism 

Nonetheless, it is relevant to explain the quantum-to-classical transition in a unified framework since 
the understanding of the decoherence (or dephasing) phenomena points out the crucial role played by the 
environmental interaction in determining whether a quantum particle shows wave behavior. Thus, there is a 
need to theoretically quantify the effect of decoherence (or dephasing) on the observed interference pattern. 
It is quite intuitive that the resulting pattern shall be an interplay between the strength of the coupling to 
the environment, the slit separation and the distance the particle travels from the slit to the screen. The 
decoherence effects on two-slit experiments have been theoretically analyzed by treating the effect of the 
environment using a phenomelogical model in [l|, [1] . In , authors described theoretically the effects on 
the interference pattern assuming the test particle develops a quantum brownian motion and solving the 
corresponding master equation, neglecting in the end the dissipation of the environment on the system. 
Contrary to these studies, authors in (lo| stated that the dynamic decoherence does not play any role in the 
visibility fringe reduction and blamed the latter on the incoherence of the source. 

In the present paper we study the visibility fringe reduction in the interference pattern of experiments with 
particles, such as fuUerenes and cold neutrons. The questions to be addressed are: how long can wc observe 
before decoherence or dephasing effects destroy the interference pattern of massive particles? How much the 
visibility fringe is reduced in these experiments and which are the possible agents of decoherence to take into 
account when modeling these experiments? Therefore, in this paper we shall study both the dephasing effects 
due to a random variable (in our case the particle's emission time) and the dynamic decoherence process 
obtained from a first principles model. Even though phenomenological models of environmental decoherence 
success fitting experimental data, we stress that a complete description of the interaction between system 
and environment is needed in order to get a well defined quantum to classical transition. 

The paper is organized as follows. In Section [H] we present the different decoherent agents that can be 
used to model this type of experiments and develop the theoretical frames to study how these agents affect 
the interference pattern. In Section IIIIl we introduce the numerical tools used in order to quantify the 
visibility fringe reduction in the pattern of an interference experiment. This is done using both analytical 
and numerical results. Section IIVI contains an application of the models described in the previous sections 
to real matter wave interferometry experiments performed with cold neutrons. Finally, in Section |Vl we 
include our final remarks. 
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II. THEORETICAL ANALYSIS 

A. Two Gaussian localized wave packets 

We shall study a typical interference experiment with particles of mass M diffracted by a grating and 
then detected on a distant screen. The particle leaves the grating and travels a distance L in the y-direction 
until it reaches the screen in a time II = ML/p^, where po is the moment's component in that direction. 
It is important to note that, in order to observe an interference pattern on the screen, particles should be 
coherent in the x-direction, whereas, the dynamics in the y-direction can be that of a free non-interacting 
particle. Hence, the experiment starts by the preparation of the initial state that emerges from the slits. 
Initially, we may reasonably assume that we have a coherent superposition of the two wave packets, centered 
at each location of the respective slits and factorized as 0, [ll[ 

^{x, 0) = (01 (x, 0) -f 02 (a^, 0)) ® xiv, 0), 

where |t/>ip and |02p correspond to the probability amplitudes for the particle to pass through slit 1 and slit 
2 (in the x-axis), respectively, while x(y,i) represents the Gaussian wave function in the y-direction (where 
no superposition is needed). Note that we are assuming translational invariance in the z-axis pj| . 

The interference pattern, in any case, corresponds to the probability distribution of the time evolved wave 
function: 

P{x,t) = (^0i(x,t)*0i(x,O + 02(x,i)*02(x,i) + 02(x,t)>i(x,t) + 0i(x,i)*02(a:,t)^|x(y,OI', (1) 

which is the diagonal part of the density matrix defined as p{x,x',t) ~ i)) t)|. 

When the system is closed, the quantum states of the system evolve accordingly to the Schrodinger 
equation. In such a case, it is easy to show that the position probability distribution on the screen at a given 
time t is: 

P{x,t)^ (^|0i(a;,t)|2-f |02(x,t)p + 2Re(0t(x,t)02(a;,t))j|x(y,t)r 

However, when the system is open, it interacts with an environment and its evolution is plagued by nonuni- 
tary features like fluctuations and dissipation, no matter how weak the coupling that prevents the system 
from being isolated is. Particularly, decoherence, as we mentioned in the preceding section, is the dynamic 
suppression of the interference terms induced on subsystems due to the interaction with an environment. For 
a superposition of locahzed wavepackets (which best describe massive particles), the initial {t = 0) four terms 
of the density matrix p{x, x' , 0) = 0i(a;, 0)*(j)i{x', 0) + (j)2{x, 0)*(j)2{x', 0)-f 02(x, O)*0i(a;', 0) + (t)i{x, 0)*(t>2{x', 0) 
correspond to four peaks. Decoherence arguments show that the off-diagonal terms die out due to the inter- 
action with the environment. As the interference pattern depends on the diagonal components of the density 
matrix, it is not obvious if the suppression of the coherences of the density matrix due to the decoherence 
process also corresponds to a disappearance of the interference pattern. In the case of open systems, the 
object of study is the reduced density matrix pr{x,x',t) of the subsystem (massive particle) which satisfies 
a master equation (sec below). 

Initially, we can assume that the environment and the total wave function of the system factorizes as 
\E'(x, 0) = [(f)i{x,0) + (f)2{x,0)]x{y,0)({X ,0), where we have introduced a new wave function C{X,t) to 
describe the state of the environment. The interference pattern at a given time t on the screen is now given 
by: 

P{x, t) = Prix, X, t)\x{y, <)P = (^\Mx, t)\^ + \Mx, t)\^ + 2T{t)Re{cl>l{x, t)Mx, t))^ \x{y, t)\^ (2) 

where T{t) encodes the information about the statistical nature of noise since it is obtained after tracing 
out the degrees of freedom of the environment. It is, in general, an exponential decaying coefficient which 
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suppresses the interference terms in a decoherence time scale . It is important to stress that in this overlap 
factor T{t) we can include not only the dynamical decoherence effects but also the dephasing ones induced 
on the subsystem due to a coupling to an external reservoir Q. 



B. Different decoherent agents 

In order to complete the analysis, we need to identify the possible "decoherent" agents so as to estimate the 
overlap factor T{t) for the different types of environment considered when modeling a two slit experiment. In 
the literature there can be found many studies that blamed the reduction of the visibility fringe on different 
causes: from the irregularities of the grating (these are named incoherence effects, and they are not really 
dynamical decoherence since they are related with the source or the preparation of the initial state) to the 
scattering of the massive particles with the air molecules to the dephasing generated by the collimation of 
the beams. 

A valid assumption, although a rather simplified version of the real problem, is the implementation of the 
model of Joos and Zeh, hereafter called scattering model in order to study the dynamics of the test 
particles moving in a quantum medium. This model, which basically is a phenomenological description of 
processes inducing loss of coherence in a quantum system, considerers that the reduced density matrix of 
the system evolves autonomously according to a markovian-typc master equation (see also Refs. jisl H^) 

dp 

i-Q^ = [H, Pr] - iA[2:, [x, pr]]. (3) 

The effect of the environment is summarized by a collision term, added to the free dynamics of the system, 
which takes into account the decoherence in the coefficient A but neglects dissipation (see discussion in 
As an example, in Rcf . , authors consider A = Aair + Aphotons, and state that the cause 
of decoherence in this type of experiments might be the scattering of the particles with air molecules and 
thermal photons during their flight from the slits to the screen . Eq. ([3]) corresponds to the many scatterer 
or high temperature limit of a more general equation [l^. Consequently, for this model, the effect of the 
environment is encoded in T(t) = exp(— At) (A is phenomenologically estimated through the wavelength 
and scattering cross section of the particles) and only considers the dynamic monito ring of the environment 
over the subsystem (i.e. dynamic decoherence base on a phenomenological model) As the contrast of 
the interference pattern is proportional to the coherence between the two paths, reduction in the contrast 
will be a direct indicator of decoherence. Therefore, spatial coherence loss of a superposition state is due to 
scattering events. In the many scatterer limit, Eq.([3]) agrees with data. Thus, decoherence is exponential 
with time and with the path separation squared, as decoherence theory usually predicts. In this context, 
other decoherence models can be applied, as the one by Hornberger, Sipe, and Arndt [T^ which uses of 
the phase space description provided by the Wigncr function to explain decoherence effects in a matter 
wave Talbot-Lau interferometer; or the more recent works by Hornberger on the formulation of the master 
equation for a quantum particle in a gas 0]. Even though we shall not considerer the Fraunhoufer limit, it 
is important to note that in Ref.[2l| thermal limitation of far-field interference has been reported. 

Another approach, which is a dephasing model, might be to consider the influence of the external classical 
time-dependent electromagnetic field on the experiment as we have previously done in . The interaction 
between the particles (electrons or neutral particles with permanent dipole moments) and classical time- 
dependent fields induces a time- varying Aharonov phase. Therein, we included a random variable to-, which 
is defined as the particle's emission time. This variable produces a fluctuating phase which averaged in 
time produces a decaying term that reduces the fringe visibility of the interference pattern. In this way, the 
uncertainty in the phase originates decoherence effects caused by the experimental difficulty of producing 
the same emission time for all particles and estimated as 

F = (e**) = lim / dtnCxp{i[Acos{u;to) + B sin{ujto)]} = M\C\), 
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where Jo is the Besscl function. The moduhis of complex number C = A + iB measures the degree of 
dephasing. The overlap factor F encodes the information about the statistical nature of noise. Therefore, 
classical or quantum noise makes F less than 1, and the idea is to quantify how slightly it destroys the particle 
interference pattern. Hence, in this case, T = F. Notably in this approach, the effect of the environment 
is constant through all the experiment since F is obtained after averaging in time and therefore, does not 
depend upon time. 

Finally, as we previously said, generally, the passage of the particles through the grating can produce 
vibrations, or other kind of interactions with the walls of the grating, able to corrupt the visibility of the 
interference pattern due to alterations in the initial coherence of the superposition. Moreover, also the finite 
size of the grating and the differences in the slit aperture can attenuate the visibility of the interference 
fringes, especially in the case of complex (very large) molecules [l9| . In the present article, we shall only 
concentrate on modeling the interaction of the interfering particles and their environment, from a microscopic 
quantum level. In this case, the dynamics of the test particles can be modeled by the quantum brownian 
motion (QBM) [2^ and the reduced density matrix of the system satisfies a master equation (see Eq.([5]) 
below) with the diffusion coefficient 'Dlt) = 2M^okBT for ohmic environment in the high temperature limit 
(when experiments at room temperature are made with large molecules, i.e. fuUerenes, and cold neutrons 
this last approximation is valid). Not only is the diffusion considered in this model environment but also 
the dissipation (through the coefficient 7(t)). Then, in this case, r(<) = cxp(— Pi) represents the noise 
induced environmental effect on the system due to the interaction with the environment. Scattering models 
or no-damped motion are just approximations obtained from our general framework [23 |. More especulative 
type of environments can be considered, such as space-time foams, quantum gravity effects, etc; but they 
are out the scope of our work since there is no experimental evidence of such decoherence agents on matter 
waves (see for example j25l.[26j). 

III. NUMERICAL ANALYSIS 

A. Interference pattern 

In this Section, we shall study the interference pattern produced by two well localized Gaussian wave 
packets, initially given by, 

^{x,0) ^ N ( cxp( ("^7 f ) + exp( + )) cM^j^-^kyy) (4) 

where 2Lq is the initial separation of the center of the wave packets, cr^g and a^^ are the initial width of 
the packet in the x and y-axis, respectively, and ky the initial moment of the particle in the y-direction. 
It is important to note that Lq, ctxO, Cyo and ky are all free parameters that have to be tuned with the 
experimental data. In addition, we assume that Apj, « py, so the moment component is sharply defined 
and the wave packet has a characteristic wavelength A^s associated A^b ~ ^/Py << ^U- 

Wc shall study the effect of decoherence on the interference pattern of an experiment with massive particles, 
by coupling our subsystem (particles) to a model environment. As we mentioned above, the experiment con- 
sists of massive particles (represented by the superposition of two localized wave packets) that are diffracted 
by a grating and registered later on a screen at a distance L. As we have already stated, the dynamics in the 
y-direction just serves to transport the particles from the slit to the screen and can then be considered as a 
"free" evolution. However, in the x-direction we need to consider a decoherent agent in order to study the 
effect of decoherence on the interference pattern observed on the screen. Thus, hereafter, we shall consider 
that the environment is a set of non-interacting harmonic oscillators and the dynamics of the test particles 
is modeled by a QBM. As noted in the preceding section, this behavior can be blamed on any interaction 
by which the particles-system become entangled with a quantum environment. In order to study the inter- 
ference pattern registered on the screen at a later time t^, we need to obtain the evolution in time of the 
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reduced density matrix pr{x, x' ,t), which is given by the following master equation 



where j{t) is the dissipative coefficient (proportional to the square of the coupling constant to the environ- 
ment), the diffusive coefficient and f{t) the coefficient responsible for the anomalous diffusion. Eq.([5]) 
has been obtained by assuming the environment to be in equilibrium, at a temperature T. In the case that 
the system is coupled to an ohmic environment in the high temperature limit {ksT >> fiio), these coefficients 
are constant ^{t) = 70, ~ 2MjokBT and f{t) « l/fcsT [2^. We restrict ourselves to the use of the 
ohmic bath since it is the type of environment which produces the correct limit for classical dissipation. It 
is the most studied case in the literature and produces a dissipative force that in the limit of the frequency 
cutoff Wcutoff — > is proportional to the velocity. In order to model more complex interactions (like charges 
with fields) it could be more appropriate to use a supraohmic spectral density. Nevertheless, it is well known 
that dynamic decoherence in the high temperature limit ocurs in a similar time-scale both for ohmic and 
supraohmic environments [23l . [27j . This is the reason why we shall only concentrate on the simplest case. 
General type of cnvironmnets can be easily included in our approach, but it is not possible to get Eq.Q as 
a limit from ^ for general non-ohmic environments. 

It is important to stress that Eq.([3]) can be obtained from Eq.(l5]) in the high temperature limit of an 
ohmic environment (neglecting dissipation) for the markovian case if written in the Lindblad form. However, 
master equation Eq. ([5]) refers to a more general movement that can be used for all temperatures and spectral 
densities, even to study the dynamics of the test particle at zero temperature (non-Markovian) limit [H, [2^. 
Yet more interesting, this formulation verifies the fluctuation-dissipation theorem for a general system in 
thermal equilibrium [24| . It is also important to note that the high temperature limit approximation is well 
defined only after a time scale of the order of l/(fcsT) 7o/2? "ensuring" the positivity of the reduced 
density matrix (a;, x', t) [ist . 

Thus, as we mentioned above, in order to study the dynamics of these two packets that best describes the 
massive particle, we need to solve the master equation Eq.([5]). The corresponding density matrix that arises 
from the initial state given by Eq.Q is 

Prix, x', 0) = Prix, x', 0) ® Priv, y' , 0) = 2N^ ( cosh(2Lo(a: + x')) + cosh(2io(a; - x'))] xiy, OTxiv', 0). 



The solution in the x-direction can be well reproduced by employing a Gaussian density matrix using the 
Born approximation [ill , [l^ . Is is worth noting that this solution does not imply far-field or Fraunhofer 
approximation. 

p,ix, x', t) = e-^(*) exp I - Ait)ix - x'f - tBit)ix^ - x'^) - C(t)(.T + x'f\ (6) 



where e"^^*-' ensures the conservation of trace. Ait) describes the range of coherence while C(i) specifies 
the extension of the ensemble in space. All functions Ait), Bit),... are real for the sake of hermicity. In our 
case, we shall study the dynamical evolution of two Gaussian wave packets located at a; = ±Lo- Therefore, 
we have to replace x x -\- and x x — L^) \n Eq. ([6]) and superpose both anzats in order to represent 
the dynamics of the two packets. In that case, the solution we shall use is 



p,ix, x',t) = 2e-^(*)e-4^«^(*) exp | - A(t)(a; - x' f - iBit)ix^ - x'^) - C(i)(a; + x')' j x 
^ cosh[4LoC(t)(.T -t- x') - i2LoBit)ix - x')] 

+ e-^^«(^(*)-^(*»cosh[4Lo^(^)(x-a;')+i2LoB(0(^ + 2;')])- (7) 
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We have numerically solved Eq.(l5]) for a free particle assuming its dynamics is modeled by QBM (in 
the X-axis) using a standard adaptative-step-size fifth-order Runge-Kutta method with initial condition 
A(0) = 1, B{0) = 0, C(0) = 1 in units of ?i = 1 = M. By doing this, we obtained the dynamic evolution 
of the coefficients A{t), Bit), C{t) and N{t). All results were found to be robust under changes in the 
parameters of the integration method. 

The intensity registered on the screen at a given time t is proportional to the position probability (diagonal 
term of the reduced density matrix) P(x,t) ^ pr{x,x,t). In the case of the initial state mentioned above, 
the intensity can be numerically obtained 

P{x, t) = e-^(*)e-4'^(*)("'-^o) cosh{8C{t)Lox) + r{t) cos{AB{t)Lox)^ , (8) 

where we have absorbed the decaying term coming from the Gaussian wave in the y-direction \xiy, i)P in the 
normalization e"^^*-*, and T{t) is a decaying exponential T{t) = e~'^^o(Mt)-Cit)) ^pj-Qvidcd A{t) — C{t) > 0). 
We know the dynamic evolution of the interference pattern at a distance L in the case the system is isolated. 
The two initial wave packets start to evolve in time and spread in the x-direction. Immediately, they start to 
develop an interference pattern. In the case of the system is interacting with a very strong environment, it is 
clear that for the same times (or even shorter ones) the interferences can not be observed because they are 
almost immediately destroyed. As the evolution continues (for a fixed length of the screen), the two packets 
continue spreading. After some time, we can no longer observe two packets on the screen because both wave 
turned into one (because of the spread of each packet). For this type of environment no interferences fringes 
will be observed for these thought experimental times. 



B. Estimation of the Decoherence Time and Fringe Visibility Reduction 



We shall estimate the decoherence time ip, i.e. the timescale for which the interferences are mostly 
destroyed, as rD(tD) = exp(— PAx^t-p) ~ 1/e. It is easily deduced that to k, l/(I?Ax^), with Ax^ = 
{x — x'Y' and T) = 2MjokBT for the ohmic environment in the high temperature limit in units of /i = 1, as 
shown on the left side of FiglTJ 
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FIG. 1: On the left: Evolution in time for the decaying exponential T{t) that destroys the interferences of the system 
for different couplings to an environment in the high temperature limit. We use units: h — c = M. = 1. Parameters 
are: Lq = 2s~^ , a^o = 0.5s~^, fcsT = 300s~^. The stronger the coupling to the environment 70 is, the sooner 
decoherence effects take place (for a fixed value of fcsT and Lq). On the right: Interference pattern registered on the 
screen at a time II for the closed and open system. In (a), we have considered the case of the isolated subsystem 
(solid red line) and the case of coupling to a strong environment with 70 = 0.01 (blue dot line). In (b), we 
have considered three different environments: 70 — 0.01 (blue dot line), 70 — 0.001 (red solid line) and 
70 = 0.0001 (black triangle line). Distance is measured in units of frequency. 
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Clearly, since the decoherence timescale depends inversely on the value of M^oksT, the stronger the 
coupling to the environment and the hotter the environment, the shorter this timescale. 

On the right side of Figll] wc can see the effects of decoherence on the interference pattern of a thought 
two-slit interference experiment with particles. In plot (a), we show the interference pattern registered on a 
screen at a distance L in a time = 0.2 s when the system is closed, i.e. there is no interaction with an 
environment, and when the system is open. In this latter case, the coupling constant is 70 = 0.01, which 
represents a strong environment because all interferences have already been destroyed (whereas they are 
present in the isolated case. In (b), we present a latter time (t^ = 0.35 s) for different coupling constants 
to the environment. We can see that for 70 = 0.01 s~^, the two wave packets are spreading and will end 
up superposing in only one final wave packet since the environment has destroyed the interference in a 
short timescale. However, for the other two environments, with smaller coupling constants, wc can see that 
the interferences are still there. Notably, the pattern remains unchanged but the visibility is considerably 
reduced as M^oksT increases. It is important to note that the visibility is considered attenuated when there 
is a lost of contrast between a maximum and a minimum with respect to the interference pattern when the 
system is isolated, i.e. the visibility is reduced when the "minimum" are not exactly zero as seen in FigU] 

At this stage, it is appropriate to quantify the loss of contrast of the interference pattern. This is done by 
defining a function called fringe visibility v, a quantity of particular importance in matter wave interferometry 

U = , 



where /max and /min represent the maximum and minimum in neighboring fringes, respectively. It is easy to 
note that the fringe visibility can be well approximated by 

^^^^ _ \pint{x,X,t)\ 



Pll{x,X,t) + p22{x,X,t) ' 

where pu = t)p, with i = 1, 2 and p-mt the interference terms. The values of this hmction range between 

(no interference fringes) and 1 (total visibility of the interference fringes). In our case, the visibility fringe 
can be numerically obtained as 



cosh{8LoC{t)x) 



Clearly, the visibility fringe goes down as t^, i.e. the observation time, is larger than the decoherence time i/j. 
However, if we succeed in performing our two slit experiment in a time < at a fixed room temperature 
ksT, we can see that the visibility fringes depends on 70 as shown on the left picture in Figl2l This is 
so, because the decoherence time depends inversely on the coupling constant. Not only can we check the 
dependence upon the coupling constant but on the separation of the slits as well. 

In the middle of Figl^l we show the visibility fringe as a function of Lq- Therein, it is clear that the 
visibility fringe goes down as the distance between the slits increases. Note that as we are plotting iy{t) for a 
fixed value of and ctxO- then wc can not vary much Lq, since wc are always assuming that a^o < Lq. We 
can also study the time evolution of the visibility i^(t), which is shown on the right side of Fig|21 Therein, we 
have plotted the evolution in time for the visibility of the first and second minimum and the first maximum of 
the interference pattern. The behavior exhibited is quite appealing. For short times, the visibility increases 
from zero to a maximum value because the interferences start to develop at that short timescale but are 
not present at t = (since the wave packets are initially separated and have to spread so as to generate the 
interferences. This maximum value coincides with the estimated decoherence time to- Then, the visibility 
starts to decrease, since the destruction of the interferences is taking place. Clearly, the decoherence is a 
dynamic process (the continuous monitoring of the environment over the test particles) and the estimated 
decoherence time is when the interferences have been reduced about a 70 %, i.e. ^{tD) ^ 1/e (see FigU]). 
However, that docs not mean that the wigncr function will be positive by that time. If one estimates the 
decoherence time as the one in which the interferences disappear completely, the estimated timescale will be 
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FIG. 2: Left: The visibility fringe v{t) plotted against 70 for a fixed II < to at room temperature (70 is in units of 
frequency). Middle: The visibility fringe i'{t) as a function of the distance between the slits Lo for a fixed time for 
ksT — 300 s~^, 70 = 0.001 s~^, Il = 0.05 s and a^co = 0.5 . All the decoherence timescales corresponding to 
the different values of Lo (in units of s^^) are checked to be longer than th- Right: Time evolution for the visibility 
fringe v{t) for ksT — 300 s~^, 70 = 0.001 s~^, Lo — 2 and Oxo ~ 0.5 . The estimation of the decoherence 
time to ~ l/(M7ofcsrLo) — 0.41 s coincides with the timescale at which the visibility starts to decrease towards a 
null value. 



longer and one might naively impute the loss of visibility on another cause but deeoherence [10|. Note that 
the visibility is a quantity that measures the loss of contrast of the interference fringes. Then, it is expected 
that those with the bigger contrast suffer from this attenuation the more, as seen in FigO Clearly, the 
observation time must be shorter than the decoherence time in order to observe the interference pattern. 
The visibility function in this case, tends to zero for longer times. 

It shall be interesting to study the visibility function for the other environmental models. In the case of 
the scattering model, the behavior of v{t) as a function of the diffusion term A and the square of the width 
of the slit Lg is qualitatively similar to that of the QBM because the expression of ^Kit) = exp(— AAx^t) 
is formally the same. Then, we expect to find that the visibility decreases as A and increases, since the 
decoherence time shall be shorter [l^l . 

Nonetheless, the visibility function for the study of the dephasing effects, i.e. when considering the in- 
teraction of the massive particle with the external time dependent electromagnetic field, is not that similar 
to the other two mentioned throughout the paper. In particular, Tc = Jo(|C|) is constant in time as we 
estimated it in (2^ for the experimental data of both neutrons and fullerenes. Therein, we calculated the 
quantity C for these massive particles and observed that, contrary to might be naively expected, in thought 
and real experiments such as the one reported in Q, Cfuiierenes ~ C'(l). However, for neutral particles with 
permanent dipole moment this value is much lower Cneutrons ^ 0(0.01) — 0(0.1). Therefore, on the left side 
of Fig[3]we present the time evolution of the visibility function lycit) defined as 

^ (t) ^ -^od^l) 

cosh{8LoC{t)x)' 

Therein, we show the time evolution of the first maximum and minimum of the interference pattern for 
different values of the C factor. It is easy to note that the development of the interferences happens in the 
same timescale of Figl2] (for the same value of Lq and axo) but in all cases, reach a different asymptotic 
value compared to the ^{t) function. The fact that fc(^) has an asymptotic limit could really be of much 
use in experiments where this effect is of importance, such as fullerenes, since once this limit is reached the 
observation time t l can be any subsequent time for the visibility function remains steady. 

Another feature of this visibility function I'cit) worth of studying is its dependence upon the separation 
of the slits Lq. On the right side of Figl3]we present this behavior. Clearly, the behavior exhibited therein 
is qualitatively different from that showed in Figj^lfor the visibility function u with T-pit). 

Finally, in Fig[4l the interference pattern for the experimental data reported in Q for two-slit experiments 
with massive particles C70 is shown. Therein, we have considered the unitary and non unitary evolution (for 
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FIG. 3: Left: Time evolution for the visibility function i^c{t) for neutrons (Cncutrons ~ 0.1) and fuUerenes (Cfuiierenes = 
1 and Cfuiieronos = 2) in the presence of an external time dependent electromagnetic field. The curves are for the first 
minimum and maximum of the interference pattern. We can see that all curves reach an asymptotic limit that is not 
zero contrary to the other two environmental models (see FigH}. Right: The visibility fringe vcit) as a function of 
the separation of the slits Lq for a fixed time for time tL = 0.20 s, th = 0.05 s and Oxo = 0.5 . The curves shown 
are for the first maximum of the interference pattern in the case of neutrons and fuUerenes. Lo is in units of . 



the three environmental models Txi{t), rA(^) and Fc) of the particles. For these massive particles, we can 
see that the interference pattern is always attenuated when the system is open. What is more significant, 
is that the effect of Fc can be as important as the other two most widely known model environments (in 
agreement with [lo| but using a different model for dephasing) and enough to model the real experiment. 
For the values of Figdl and asking to > ti, (and correspondingly > t^), we obtain a constraint for the 
free parameters of each model: 70 < 7.14 x 10~^ [s~^] (as estimated in and A < 7.44 x 10^^ [m~^s~^] 
(approximately the value used in [ll|) for the experimental data at room temperature reported in 

We want to emphasize that Eq.© phenomenologically models the decoherence effects neglecting the 
dissipative process. The value we obtained for 70 is extremely small so a valid question might be if it is 
necessary to include dissipation in the model. We state that it positively is in order to have a complete 
and formally correct description of the process. By including the term proportional to 7(t) in Eq.(l5]) we are 
assuring the fulfilment of the fluctuation-dissipation theorem, also known as Einstein formula in the high 
temperature limit. It is known that the diffusive coefficient V = 2M^okBT is proportional to 70. In this way, 
if 7o happens to be zero (which means no dissipation), the diffusive term would also be zero. The correctness 
of the formulation can also be checked in the fact that even though 70 is extremely small, M^^kbT can be very 
large. In such a case, the decoherence effects would be very important whereas the dissipative interaction 
between the particles and the environment can be ignored. In other words, small dissipation implies that 
the particles could have a neglible damping term in the semiclassical Langevin equation of motion along the 
X direction, but the existance of noise ensures that decoherence shall be effective. 

It is important to stress that all these environmental models consider one and only one "decoherent" agent 
influencing the interference experiment. However, all these effects can be together considered to be present 
in a two-slit experiment. In such a case, the attenuation factor F would be F w Fa + Fp -|- Fc. The effect 
of the environment would be equal to the sum of the three factors (as seen in Figl4]). Therefore, it is enough 
to consider the biggest one (unless they are all the same order of magnitude). All in all, it is necessary to 
remark that decoherence effects do play a crucial role in the fringe visibility reduction. On the one hand, 
authors in said that the incoherence of the source was to blame for the fringe visibility reduction. On 
the other hand, authors in developed a phenomenological theoretical model where decoherence was a 
priori introduced by assuming an exponential damping of the interferences. So far, we have shown that 
the scattering of the massive particles with the air molecules and dephasing (for example introduced by the 
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FIG. 4: The interference pattern {v ^ 0.68) registered on the screen for the unitary evolution and the nonunitary 
evolution considering the different model environments for a double slit diffraction experiment with massive particles 
C70. The curves are done with the experimental data reported in 0]. The values used for the plot: 70 = 2.6x 10~* s~^, 



random emission time) of the experimental setup are all responsible for the fringe visibility reduction and 
approximately of the same order of magnitude. 



IV. APPLICATION: EXPERIMENTAL DATA FOR NEUTRONS 

In this Section, we shall use the existing experimental data [sol . [sH to reproduce the observed patterns for 
neutral cold atoms. 

If we consider that II = L/v ~ MX^^L/ [2T:h) and » MaxoLo/h, then the position distribution on 
the screen at this time II can be well approximated by: 

where ^{^l) depends on the model environment we want to use to describe the conditions in which the two- 
slit experiment is being done evaluated in the observation time. In this way, Eq.Q describes the intensity on 
the screen as a function of the experimental parameters, i.e, the mass M of the cold neutrons, the associated 
wavelength AdB, the distance to the screen L, the distance between the slits (assuming the two slits are as 
similar as possible) and the initial width of the wave packet Uxo- AH these values can be found, for example, 
in Q for cold neutrons. Note that Eq.Q is equivalent to Eq.® , identifying our time dependent theoretical 
parameters with the real experimental ones. Thus, we have for a fixed observational time th (making the 
same assumptions as in the above section) 

= — F ^"^"^ = — T~ ■ 

In the case we studied in the preceeding Section, assuming that the dynamics of the test particle can 
be modeled by a quantum brownian motion, the expression for r(ii) is Vj:,(ti^) = cxp^—tL/to) with to = 
VIT? j (AI^okBTLQ) where we have reincorporated h. In the estimation of this time we have considered that 
Ax^ ~ Lg, which in fact is an underestimation of the decoherence time for lengths bigger than Lq. As 
the experiment is done at room temperature, the only free parameter is the value of 70. On the other 
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FIG. 5: The interference pattern {v ^ 0.57) registered on the screen for the unitary evolution and the nonunitary 
evolution considering the different model environments. Experimental results obtained by Zeilinger et. al [s^l for 
the double slit diffraction of cold neutrons are plotted with blue star dots. The values used for the plot: 70 = 
5.0 X IQ-i^ A = 5.5 X lO" s'V 

hand, if the model environment were assumed to be the one of the scattering with air molecules, where 
the effects of the environment are included in the collision term A, then the expression for r(ti) would be 
T^itL) = exp(-tL/tA) with tA = 3/(AL2). 

The other possible environmental model mentioned in Scc lIIBI was to considered the interaction with the 
charged or neutral (with permanent dipole moment) particles with the electromagnetic field. This is an 
interaction that is always present, can never be turned off, although sometimes it is possible to neglect it. 
As we previously studied in [2^ . in the case of neutral particles with permanent dipole moment, this effect 
is not so important. 

In Figl5]we have plotted, the interference pattern for the experimental data reported for experiments with 
cold neutrons for the isolated and open system. In this last case, the nonunitary evolution has been considered 
for the different environment models mentioned above. We can clearly sec, that for this case, the nonunitary 
evolution when considering the interaction of the cold neutrons with a time-varying electromagnetic field 
(orange dots) is exactly superposed with the unitary one (dotted line). That means that the incoherence 
effects can be completely neglected. However, the other two model environments, whose effects are considered 
in Tf^{t) (black triangle dotted line) and Vpit) (red solid line), fit correctly the experimental data obtained 
by Zeilinger et al. (blue star dotted line) in [32| . The fact that the interference pattern is observed implies 
that the decoherence time is Id (and tjs) is larger than the observation time II- That sets us a constraint 
to the expected values for 70 (and A), the free parameter in each model. By asking to > t^, we have 
7o < 8 X 10~^ In thfi case of modeling the environment by a collision term A, if t\ < ti, is asked, then 

A < 1.28 X IQ^^ [m-2s-i]. 

V. FINAL REMARKS 

The effect of the environment on the interference pattern of a two-slit interference experiment with massive 
particles has been studied phenomenologically in the literature. 

However, here we have presented a fully quantum mechanical treatment using a microscopic model of 
environment and also a concrete example to include dephasing effects. Therefore, we have studied the effects 
of decoherence on the interference pattern of thought experiments and presented an analysis of matter wave 
interferometry in the presence of a dynamic quantum environment such as the quantum brownian motion 
model. We have shown the interference patterns and visibility function v[t) for thought diffracted free 



13 



particles and analyzed their dependence upon different parameters of the model in the high temperature limit 
(assumption valid for massive particles interfering at room temperature). As was expected, the visibility 
decreases as the value of the diffusion coefficient increases and, in particular, we showed this effect for 
different values of the coupling constant 70 to the environment. What is more important, we have seen 
that the visibility fringe is considerably reduced when considering an open quantum system, although the 
structure of the interference pattern remains unchanged. 

Yet more important, we defined the visibility function vc for a model environment previously developed 
which describes dephasing effects originated in the experimental difficulty of producing the same initial/final 
state for all particles (i.e the existence of a random variable such as the particle's emission time). We showed 
that it is qualitatively different than the one commonly found in the literature and very important in the 
case of experiments with massive particles such as fullerenes. In this case, dephasing effects are enough to 
model the attenuation of the interference pattern observed in the real experiment, whereas in the case of 
cold neutrons these effects are not of such importance. Therefore, in the latter case we must consider the 
decohcrence effects by using the corresponding formulation. This result might have been expected since the 
interaction of more massive particles with the external classical field is more important than for those with 
a smaller mass where other kind of interactions seem to prevail. 

Finally, the effect of the environment on a two-slit experiment can be modeled by considering different 
effects such as the scattering of the massive particles with the air molecules, the randomness of the arrival 
or emission times and the presence of a classical time dependent electromagnetic field. Even though there 
exist conceptual differences in all the cases mentioned throughout the paper, we showed that all these effects 
reduce the visibility fringe and can be formally deduced from a microscopic model (whether the QBM for 
decoherence effects studied in this paper or a fluctuating Aharonov-Casher phase studied in our previous 
contribution) . They are all included in the noise induced effects introduced in the subsystem when the latter 
is coupled to a quantum external environment. 
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